ON PROJECTIVE SYSTEMS OF RATIONAL DIFFERENCE 

EQUATIONS 

FRANK J. PALLADINO 

Abstract. We discuss first order systems of rational difference equations which have 
the property that lines through the origin are mapped into lines through the origin. 
We call such systems projective systems of rational difference equations and we show 
a useful change of variables which helps us to understand the behavior in these cases. 
We include several examples to demonstrate the utility of this change of variables. 



1. Introduction 

We study first order systems of rational difference equations which have the property 
that lines through the origin are mapped into lines through the origin. We call such 
systems projective systems of rational difference equations. Consider the most general 
first order system of linear fractional rational difference equations. 

a l 4" X/i=l 01,i x n,i t\j 

Xn + i,i = k , n G N, 

"2 + Ei=l faiX n ,i _ TXT 

xn + i,2 = k , n G N, 

A 2 + 2^i=l B 2,i%n,i 



x n +i,k = I , n e N. 

■^■k / jj—i -L>k,i%n,i 

It is customary to assume that initial conditions and parameters are nonnegative. We 
will also make this assumption. However, we note that if proper care is taken to avoid 
problems which may arise from divison by zero, then the changes of variables we present 
in this article will also work in the case of complex parameters and complex initial 
conditions. In the following three situations this system is projective. If on = = A{ 
for all 1 < i < k then the system is projective. If ol\ = for all 1 < i < k, — Aj 
for all i,j G {1, . . . , k}, and = Bij for all G {1, . . . , k}, then the system is 

projective. If A, = for all 1 < i < k, cti = aj for all i,j G {1, . . . , A;}, and fin = fiej for 
all I G {1, . . . , k}, then the system is projective. 
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In words, the three scenarios can be described as follows. The first scenario is when 
all of the constants in the numerator and denominator, cnj and A*, are zero. Such sys- 
tems of rational difference equations are commonly referred to as homogeneous systems 
of rational difference equations. Homogeneous systems of two rational difference equa- 
tions have been studied extensively in [3]. The strategy employed in [3] was to use the 
projective change of variables to carry out a reduction in the number of systems, from 
two to one. This allowed the authors of [3j to determine the behavior in many cases of 
rational systems in the plane. We will call such systems projective rational systems of 
homogeneous type. 

The second scenario is when all of the constants in the numerator, are zero and 
the denominators are all the same. In this scenario, the projective change of variables is 
commonly used to embed the linear fractional system into a higher dimensional linear 
system. An example of this is the solution of the Riccati difference equation which can 
be found in [2] and pE]. We will call such systems projective rational systems of linear 
type. 

The third scenario is when all of the constants in the denominator, Aj, are zero and 
the numerators are all the same. We will call such systems projective rational systems 
of hyperbolic type. We will focus on this scenario in several examples to show how these 
changes of variables can be used to simplify the problem considerably and determine 
global behavior. 

In the next section of this article, we provide a change of variables for each of the 
three types of projective rational difference equation. We then, in the following section, 
give examples for each type of projective rational difference equation. We show how the 
change of variables can be used to determine the qualitative behavior of rational systems 
in several special cases. 



The general projective rational difference equation of homogeneous type is the follow- 
ing system of rational difference equations. We will assume positive initial conditions 
for our work here. 



Now we divide through by one of the variables, in this example we will use x n+ i^ but 
this choice is arbitrary and any varaible may be used. In order to do this, we should 
have k > 1. We get 



2. A CHANGE OF VARIABLES FOR EACH SCENARIO 





ti=l Pk,i%n,i 

k Bi. x ■ 



n E N. 
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Ek a x n ,i \ I v^fc tj x Uti \ 

■■„T,..-, M"*-"^ ^1%^ 

(EL 1 A,.S)(El 1 B.-i,ti)' 



Using the change of variables u n j = we get the following. 



(Pl,k + Ei=l Pl,i u n,i) \ Bk,k + Ei=l -^fc.iMri 
Mn+1,1 = y — y, — r-' H E N ' 



Un+l,fc-l = -7 — yj —± V' nEN - 

— \,iU"n,i ) 

The general projective rational difference equation of linear type is the following system 
of rational difference equations. We will assume positive initial conditions for our work 
here. 

s ra +ii= S^lfa^! n G N, 

^1 + Li=l 5 l,^n,i 



A l + Ei=l B^Xnj 

Now we divide through by one of the variables, in this example we will use x n+ i^ but 
this choice is arbitrary and any varaible may be used. In order to do this, we should 
have k > 1. We get 

Ek n Xn,i 
i 1 ' -* ! /i.A- 



nGN, 



re ge 

Pfc,fc + Z^i=l Pk,iV n ,i 



Using the change of variables v n j = we get the following. 
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Pk-l,k + Yli=l (3k-l,iV ni 

v n+ i,k-i = — ; ^ fc _i ~ , n E N. 

Pk,k + Z^i=l Pk,iVn,i 

The general projective rational difference equation of hyperbolic type is the following 
system of rational difference equations. We will assume positive initial conditions for 
our work here. 

%n+i i = 1 — — — ' ~, n G N, 

V Bi x ■ 



£n+i fc = i — — — : — n E N. 

/ yj— ] 1 -'k,i Jj n,i 

Now we divide through by one of the variables, in this example we will use £ n +i,fc but 
this choice is arbitrary and any varaible may be used. In order to do this, we should 
have k > 1. We get 



.1 



n+1,1 



n e N, 



n E N. 



Using the change of variables w n j = we get the following. 

B k ,k + J2i=i B k ,iW ni 
w n+i,i = ~ ^s^k-ij, > n EN, 

-Dl,k + Aui=% J=>l,iW n ,i 



Bk,k + Yli=i Bk,iiv n i _ 



Bk-i,k + Si=i Bk~i,iW nt i 

3. Examples and remarks 

We begin with a well known example which demonstrates the projective change of 
variables of linear type. Our first example provides a different change of variables for 
the Riccati equation than the standard change of variables found in [I]. 

Example 1. Consider the difference equation 

a + (3x n _ 
A + Bx n 
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with initial condition Xq G (0,oo) and parameters a,fi,A,B G [0, oo) such that (A + 
B)(a + ft) 7^ 0. Inspired by the projective change of variables of linear type we may 
construct the following linear system in the plane, 

y n+1 = az n + (3y n , n = 0,1,2,..., 

z n+1 = Az n + By n , n = 0, 1, 2, . . . . 

Choosing initial conditions y$ = Xq and Zq = 1 we have, — = x n for all n > 0. We prove 

Zn 

this by induction, n = provides the base case. Notice that since (A + B)(a + pi) ^ 
and 2/o, %a G (0, oo), y n , z n G (0, oo) for all n G N. Now, suppose j 1 = x n , then 

« + fan _ 01 + _ ®Z n + ftyn _ Vn+1 

Xn+l ~ A + Bx n ~ A + B^~ Az n + By n ~ z n+l ' 
So — = x n for all n > 0. 

In our next example, we use the projective change of variables of hyperbolic type 
and the results from [5] to determine the qualitative behavior of a system of 3 rational 
difference equations. 

Example 2. Consider the following system of three rational difference equations, 

%n+l — ~T~\ ' ^ — 0, 1, 2, ... , 

Vn+1 = 7^—, 71= 0,1,2,..., 

Dz n 

Zn+1 = T, f ' n = M,2, 

px n -\- OiZ n 

with initial conditions xo, yo, zq G (0, oo) and parameters C, A, D, ft, a G (0, oo). For this 
system, {x n }^L 1 converges to 

2 (/3 D -AD- a) 2 + AaDf3C + 2 ({3D — AD — a) a/ (f3D -AD- a) 2 + 4aD/3C 



x 



(2/3C) (2af3C + f3 2 D - AD[3 - a/3 + Py/(/3D -AD- a) 2 + 4aD/3<7) 

{y n }^ =1 converges to 

__(3D-AD-a+ y/(pD -AD- a) 2 + 4aDf3C 
V ~ 2(3CD ' 

{z n }™ =1 converges to 

/3D-AD-a + y/(pD -AD- a) 2 + 4aD(3C 
2af3C + f3 2 D - AD{3 - af3 + fty/IfiD -AD- a) 2 + 4aDf3C' 
Proof. Dividing through by z n+ i we have the following relationship, 

x n+ i Px n + az n + a 



Zn+i Cy n + Az n C^ + A'' 



n = 0,1,2, 



Vn+1 P x n + az n [3x n a 

= = 1 , ra = 0,l,2, 

Z n +l 

Dz n Dz n D 
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Relabeling so that u n = — and v n = — we have, 

(3u n + a 

u n +i = y; —7, n = 0,1,2,..., 

Cv n + A 

(3u n a 
v n+1 = —- + —, n = 0,1,2,.... 

So the sequence {u n }^ =l satisfies the following second order rational difference equation, 

(3u n + a 



~ ffC Mn - 1 . . . r, ' ™ — J-j^} • • • ) 



with initial conditions iio = and Ui = We know from [5] that {ti n }^i converges 

to 



u 



[3D — AD — a + V /(^D - AD - a) 2 + 4aD/3C 



2(3C 

Thus, {v n }'%L l converges to 



2aC + (3D - AD - a + J ((3D -AD- a) 2 + AaDfiC 

v = . 

2CD 

Thus, since y n+1 = ^, {y n }n=i converges to 



x = 



__(3D-AD-a + a/ ((3D -AD- a) 2 + AaD(3C 
V ~ 2f3CD ' 

Since x n = ^ IL , {x n }^ ( L 1 converges to 

2 ((3D -AD- a) 2 + AaDpC + 2 ({3D — AD — a) ^(fiD -AD- a) 2 + 4aD(3C 
(2/3C) (2af3C + f3 2 D - AD(3 - a/3 + (3^/{(3D -AD- a) 2 + 4aD/3<7) 

Since z n = — , \z n }°°_, converges to 

(3D - AD - a + y/((3D -AD- a) 2 + 4aDf3C 
Z ~ 2a(3C + (3 2 D - AD(3 - a/3 + /3y/(/3D -AD- a) 2 + 4aD(3C 

This example demonstrates the utility of the projective change of variables. If we 
were to try to analyze this system of three rational difference equations without using 
the projective change of variables and Merino's results in [5], then we would be forced 
to overcome the same obstacles already overcome by Merino with the help of a lengthy 
Mathematica computation. □ 

In our next example, we use the projective change of variables of homogeneous type 
to reduce a 3 dimensional system which is neither competitive nor cooperative, to a 
cooperative system in the plane. We then use the available theory to determine the 
behavior for the reduced system. This allows us to determine the qualitative behavior 
of the original system of 3 rational difference equations. 
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Example 3. Consider the following system of three rational difference equations, 

_ %n + Vn _ n 1 9 

x n+l — a , 71 — U, ±, Z, . . . , 

■fi-iZn ~\~ x n -\- y n 
Vn+i = A I — — - , n = 0,1,2,..., 

^n+i = , n = 0,1,2,..., 

with initial conditions xo,yo, z G (0, oo) and positive parameters a, A\, and A 2 . Define 
a polynomial, 

P(w) = -w 3 + (2 - Ai - A 2 - a) w 2 + (A x + A 2 - o^4i - aA 2 - A X A 2 ) w - aA 1 A 2 . 
Further, define 



2(2-A 1 -A 2 -a) + a/4 {2-A 1 -A 2 - af + 12 (A t + A 2 - aA 1 - aA 2 - A 1 A 2 ) 



w r , 



Whenever P{w rn ) < or one of the following two conditions hold, 

(1) Ax + A 2 - aAi - aA 2 - A X A 2 < and 

2 (2 - A x - A 2 - a) < y/-12 (A~i + A 2 - aA x - aA 2 - A X A 2 ), 

(2) Ai + A 2 - aA x - aA 2 - A X A 2 = and A 1 + A 2 + a > 2, 
then 

lim n ^ QQ x n 0, 

li^n^ooUn = 0, 

Whenever P{w m ) > and one of the following three conditions hold, 

(1) Ai + A 2 - aAi - aA 2 - A X A 2 < and 

2 (2 - A ± - A 2 - a) > y/-12 (A t + A 2 - uA x - aA 2 - A t A 2 ), 

(2) A + A 2 - aAi - aA 2 - A X A 2 = and A x + A 2 + a < 2, 

(3) A ± + A 2 - aA x - aA 2 - A X A 2 > 0, 

then the equation 

w 2 ( 1 1 

W = = ~A = + 'A 

a + w \A X + w A 2 + w / 
has three nonnegative solutions w — 0,w 1 ,w 2 with < w 1 < w 2 . 
If Bi±m e[0}Wlh then 

li7n n ^. OCl x n 0, 
lim n -^ 00 y n = 0, 

7/5a±aa =Wl} then 

- Wl •> i 

Ai + wi 



7/5a±aa e (^,00), tfien 
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yn = , , _ , n>l, 
A 2 + t«i 

a 

z n = 1 + — , n > 1. 



lim n ^ ! . 00 y ri 



Ax + w 2 
w 2 

A 2 + w 2 

_ a, 

w 2 

Proof. Dividing through by z n+ i we have the following relationship, 

2 

X n I Vn 



X n +1 



-, n = 0,l,2,..., 



+i ( Ai + ^ + mA (a+^+ V -A 

, n = 0,l,2,.... 



In I Jn 

Vn+1 



z n+l 



Relabeling so that u n = — and v n = — we have, 



%L n -\-\ = ~r~: r~~T r, 77/ = U, 

(Al + U n + t> n ) (a + U n + V„) 

(ti n + f n ) 2 

u n+1 = — — -, n = 0,l,2, 

(A 2 + u n + v n ) (a + u n + w„) 



w 2 , / 1 1 



Let w n = u n + v n , then 

w n+ i = — ^- ( - - + - - ) ■ // = 0.1.2.... 

a + w n \Ai+w n A 2 + w n 

This is a first order difference equation in one variable where the function is monotone 
increasing, so we may use the stair step diagram to find the qualitative behavior of 
solutions. First, let us find the equilibria. We solve 

w 2 ( 1 1 



w 



a + w \Ax + w A 2 + w J 



_ f w (Ax + A 2 + 2w) -(a + w) (A 1 + w) (A 2 + w) 

~ \ ^ a + w )(A 1 + w) (A 2 + w) 

q_ w [ ~ w3 + (2 - A - A 2 - a)w 2 + (Ax + A 2 - aA x - aA 2 - A X A 2 ) w - aA x M 

~ W \ ( a + w) (A x + w) (A 2 + w) 

So w = is an equilibrium and by Descartes' rule of signs there are either 2 positive 
equilibira, or zero positive equilibria. Let us examine the polynomial, 

P(w) = -w 3 + (2 - A x - A 2 - a) w 2 + (Ax + A 2 - aAx - aA 2 - AxA 2 ) w - aAxA 2 , 
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more closely. The derivative of this polynomial is 

D(w) = -3w 2 + 2 (2 - A x - A 2 - a) w + A x + A 2 - aAx - aA 2 - A X A 2 . 
So, if A-y + A 2 - olA x - aA 2 - A X A 2 > 0, then the root 

2(2-A 1 -A 2 -a) + J A (2 - A x - A 2 - of + 12 {Ax + A 2 - aAx - aA 2 - AxA 2 ) 



is the only positive root of D(w). Since D(0) > and Urriyj^ooDiw) = — oo in this case, 
w m is a local maximum for P(w) and the largest value of P(w) for positive w. Since 
P(0) < 0, we have two possibilities in the case where Ai + A 2 — aAi — aA 2 — A\A 2 > 0. 
If A\ + A 2 — aA\ — aA 2 — A\A 2 > and P(w m ) > 0, then there are three equilibria 
w = 0,Wi,W 2 with < W\ < w 2 . In this case, 



is negative on the intervals (0,Wi) and (lo 2 ,oo) and positive on the interval (Wi,w 2 ). 
So, we have that the basin of attraction of is [0,Wi), the basin of attraction of W\ is 
{Wi}, and the basin of attraction of w 2 is (Wi, oo). 

If Ai + A 2 — aA\ — aA 2 — A\A 2 > and P(w m ) = 0, then this contradicts Descartes' 
rule of signs since there is a single positive equilibrium in this case. If A\ + A 2 — olA\ — 
aA 2 — AiA 2 > and P(w m ) < 0, then the function 



is negative on (0, oo) so there is only one non-negative equilibrium w = and this 
equilibrium is globally asymptotically stable. 

If A\ + A 2 — aAi — aA 2 — A\A 2 = and A\ + A 2 + a > 2, then by Descartes' rule 
of signs P(w) has no positive zeros, so f{w) is negative on (0, oo) so there is only one 
non-negative equilibrium w = and this equilibrium is globally asymptotically stable. 

Now, whenever Ai + A 2 — aA 1 — aA 2 — A\A 2 = and A\ + A 2 + a < 2, then 



and lim w ^ > . 00 D{w) = — oo. Thus, in this case, w m is a local maximum for P{w) and the 
largest value of P(w) for positive w. Since P(0) < 0, we have two possibilities in the 
case where A x + A 2 - aA x - aA 2 - A X A 2 = and A ± + A 2 + a < 2. UA 1 +A 2 - aAi ~ 
aA 2 — A\A 2 = 0, Ai + A 2 + a < 2, and P(w m ) > 0, then there are three equilibria 
w = 0,Wi,w 2 with < Wi < w 2 . In this case, f(w) is negative on the intervals (0,Wi) 
and (uJ2,oo) and positive on the interval (Wi,w 2 ). So, the basin of attraction of is 
[0, Wi), the basin of attraction of Wi is {wj, and the basin of attraction of w 2 is (Wi, oo). 

If Ax + A 2 - aAx - aA 2 - A X A 2 = 0, Ax + A 2 + a < 2, and P(w m ) = 0, then this 
contradicts Descartes' rule of signs since there is a single positive equilibrium in this 
case. If Ax + A 2 - aAx - aA 2 - A X A 2 = 0, Ax + A 2 + a < 2, and P(w rn ) < 0, then the 
function f(w) is negative on (0, oo) so there is only one non-negative equilibrium W = 
and this equilibrium is globally asymptotically stable. 



w. 



3 
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If A 1 + A 2 - olA x - aA 2 - A X A 2 < and 



2 (2 - A x - A 2 - a) < y/-12 (Ax + A 2 - aAx - aA 2 - A X A 2 ), 

then D(w) has no positive real roots. Since D(0) is clearly negative, we get that D(w) 
is negative on (0, oo). Thus, since P(0) is negative, P(w) is negative on (0, oo). Thus, 
f(w) is negative on (0, oo). So there is only one non-negative equilibrium w = and 
this equilibrium is globally asymptotically stable. 
If Ai + A 2 - aAi - aA 2 - A X A 2 < and 

2{2-A 1 -A 2 -a)= y/-l2 {Ay + A 2 - aA x - aA 2 - A 1 A 2 ), 

then D(w) has only one positive real root. Since D(0) is clearly negative and lim w ^oo = 
— oo, we get that D(w) is nonpositive on (0,oo). Thus, since P(0) is negative, P(w) is 
negative on (0, oo). Thus, f(w) is negative on (0, oo). So there is only one non-negative 
equilibrium w — and this equilibrium is globally asymptotically stable. 
If A x + A 2 - aA x - aA 2 - A X A 2 < and 



2 (2 - Ax - A 2 - a) > y/-12 (Ax + A 2 - aAx - aA 2 - AxA 2 ) } (1) 
then D(w) has two positive real roots. D(0) is clearly negative and 
^ . 2 - Ax - A 2 - a\ (2 - A x - A 2 - a) 2 

D = ^ '- -Ax + A 2 - aAx - aA 2 - A X A 2 . 



3 

Thus by the inequality 



Moreover, lim^^oo = — oo. Thus, in this case, w m is a local maximum for P(w) and 
either P(0) or P(w m ) is the largest value of P(w) for nonnegative w. Since P(0) < 0, we 
have two possibilities in the case where A x + A 2 — aAx — aA 2 — AxA 2 < and (1) holds. 
If Ax + A 2 — aAx — aA 2 — A X A 2 < 0, (1) holds, and P(w m ) > 0, then there are three 
equilibria w = 0,Wx,w 2 with < w x < w 2 . In this case, f(w) is negative on the intervals 
(0,Wx) and (w 2 , oo) and positive on the interval (Wx,w 2 ). So, the basin of attraction of 
is [0, wTi), the basin of attraction of Wx is and the basin of attraction of w 2 is 

(wx, oo). 

If Ax + A 2 — aAx — aA 2 — AxA 2 < 0, (1) holds, and P(w m ) = 0, then this contradicts 
Descartes' rule of signs since there is a single positive equilibrium in this case. If Ax + 
A 2 — aAx — aA 2 —AxA 2 < 0, (1) holds, and P(w m ) < 0, then the function f(w) is negative 
on (0, oo) so there is only one non-negative equilibrium w = and this equilibrium is 
globally asymptotically stable. 

Thus summarizing the results obtained above via the stair step method from [4], we 
get the following result. Whenever P(w m ) > and one of the following three conditions 
hold, 

(1) Ax + A 2 - aAx - aA 2 - A X A 2 < and 



2 (2 - Ax - A 2 - a) > y/-12 (Ax + A 2 - aAx - aA 2 - A X A 2 ), 
(2) Ax + A 2 - aAx - aA 2 - A X A 2 = and Ax + A 2 + a < 2, 
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(3) A x + A 2 - aA x - aA 2 - A X A 2 > 0, 
then there are three equilibria w = 0,Wx,w 2 with < W\ < w 2 . The basin of attraction 
of is [0,Wx), the basin of attraction of w\ is {w{\, and the basin of attraction of w 2 is 
(Wi, oo). 

Since w n = Xn+y ™ we have: 

x n+i — Al j n — 0,1,2, ... , 

1 

Vn+l = M +l ' n = 0,1,2,..., 

Wn 

a 

z n+1 = H , n = 0, 1,2, . . . . 

Thus in this case if 2^±Mi e [o,^), then 

liTn n — yoc%n 

l'Wi n ^ QC :y n = 0, 

/ 2 Tfln — v 00 'S-n OO. 



If 2^ = «Ji, then 



If ^±m e (iy 1>0 o), then 



- Wl > 1 

2/n = , — , n > 1, 
A 2 + t«i 

a 

z n = 1 + — , n > 1. 



w 2 



A 2 + «y 2 ' 

U>2 



Whenever P(w m ) < or one of the following two conditions hold, 
(1) A x + A 2 - aAi - aA 2 - AiA 2 < and 



2 (2 - A ± - A 2 - a) < y/-12 (Ax + A 2 - aA 1 - aA 2 - AxA 2 ), 

(2) Ax + A 2 - aAx - aA 2 - A X A 2 = and A 1 + A 2 + a > 2, 
then 

liui n —^ 00 x n 0, 
liui n ^. 00 y n = 0, 



□ 
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In our next example, we use the projective change of variables of hyperbolic type to 
reduce a 3 dimensional system to the rational system in the plane numbered (11,11), 
in the numbering system introduced in [T]. We then further reduce the system to a 
second order difference equation which decouples into two Riccati difference equations 
for which the solutions are known. This allows us to determine the qualitative behavior 
of the original system of 3 rational difference equations. 

Example 4. Consider the following system of three rational difference equations, 

x n+i — ; > Ti — 0, 1, 2, . . . , 

Az n + By n 

yn+i = Cz T Dr > n = 0,l,2,..., 

1 

z n +i = — , n = 0, 1,2, 

Zn 

with initial conditions Xo,yo,Zo G (0, oo) and parameters A,B,C,D G (0,oo). Then 
every solution converges to the not necessarily prime period 2 solution, 



X2n — Zq 



U2n = Zq 



X2n+1 



V2U+1 



D - 


AC - 


B 


+ ^{D-AC- 


B) 2 + AADC 








2AD 










2A 




AC 


+ D- 


B 


+ y/(D -AC- 


B) 2 + AADC' 








Z2n — Zq, 




D — 


AC- 


B 


+ y/(D-AC- 


B) 2 + AADC 








2ADz 










2A 





1 

Z2n+1 — — • 
ZQ 

Proof. Dividing through by z n+ i we have the following relationship, 

x n+l 1 



z n+l A + B^- 



, 71 = 0,1,2,..., 



Vn+l 1 n 1 o 

Relabeling so that u n = — and v n = — we have, 

1 

= a I p ' n = 0,1,2,..., 
A + Bv n 
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V n+1 = ——^— , 72 = 0,1,2,.... 

C + Du n 

So the sequence {it n }^i satisfies the following second order rational difference equation, 

with initial conditions uq = u and U\ = A + Bv ■ This equation decouples into two Riccati 
difference equations. Thus every solution converges to the positive equilibirum 

__D-AC-B + y/(D — AC — B) 2 + AADC 
U ~ 2AD • 

Thus {v n }^ =1 converges to the positive equilibirum 

2A 



AC + D - B + y/(D - AC - B) 2 + AADC 

Now, clearly z 2 n = z and z 2n +i = ^ for the sequence {z n }^ =1 . Since converges 
to the positive equilibirum 



u = 



D-AC-B+ y/(D — AC — B) 2 + AADC 



2AD 

and u n = J 1 , {x n }^i converges to the period two sequence, 

D-AC-B + y/(D -AC- B) 2 + AADC 

and 

D-AC-B+ ^{D - AC - B) 2 + AADC 
X2n+1 ~ 2ADz~ • 

Since {v n }™ =1 converges to the positive equilibirum 

2A 



AC + D-B+ y/(D - AC - B) 2 + AADC' 
and v n = I 2 -, {yn}^Li converges to the period two sequence, 

2A 



V2n Z °AC + D-B+ ^{D -AC- B) 2 + AADC' 

and 

2A 

V2n+l — 



(AC + D-B+ y/(D -AC- B) 2 + AADC^j 



z 



Thus, for our original system, every solution with initial conditions x ,y ,z , converges 
to the not necessarily prime period 2 solution, 



D-AC-B + y/(D -AC- B) 2 + AADC 

*2n ~ Z — , 
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2A 



V2n = Z 



AC + D 



B + y/(D -AC- B) 2 + AAD& 



Z2n — Zq, 



%2n+l 



D-AC 



B + y/(D -AC- B) 2 + AADC 



2ADz 



2A 



V2n+1 



(AC + D-B+ y/(D -AC- B) 2 + AADC^j ' 



1 



Z2n+l 



□ 



4. Conclusion 



We have presented 3 full families of projective rational difference equations with the 
corresponding changes of variables. This should be of interest to researchers in the area 
for two reasons. First, we point out the geometric origins of well known changes of 
variables. Second, we present new families of projective rational difference equations 
and give a useful change of variables in each case. 

Moreover, we have presented several special cases where the projective change of 
variables can be used to completely determine the qualitative behavior. These special 
cases are interesting in their own right and it would have been exceedingly difficult to 
determine the qualitative behavior in certain cases without the projective change of 
variables. 
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